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Discovery of a 30 day period in the supergiant fast X-ray transient 
SAX J1818.6-1703 
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ABSTRACT 

SAX J1818. 6-1703 has been characterised as a Supergiant Fast X-ray Transient sys- 
tem on the basis of several INTEGRAL/lBlS detections since the original BeppoSAX Wide 
Field Camera detection. Using IBIS/ISGRI, Swift/BAT and archival observations, we show 
that in fact SAX J1818. 6-1703 exhibits emission on a period of 30 ± 0.1 days, with a high 
degree of recurrence. SAX J1818. 6-1703 is therefore the second SFXT shown to exhibit pe- 
riodic outbursts, but with a considerably shorter period than the other known system, IGR 
J11215-5952. 

Keywords: gamma-rays: observations - X-rays: binaries - X-rays: individual: SAX J 18 18.6- 
1703 



1 INTRODUCTION 



In just over 5 years, the IBIS instrument dUbertini e t al."2003') on 
board the INTEGRAL gamma-ray satellite jwinider et al. 2 0031) 
has revolutionised our classical view of High Mass X-ray bina- 
ries (HMXB). Specifically, IBIS has discovered many new HMXBs 
hosting massive OB supergiant stars (the so-called Supergiant High 
Mass X-ray binaries, SGXB); in just a few years their population 
has been almost tripled thanks to the IBIS identification of two pre- 
viously unrecognized su bclasses of SGXBs: the highly obscured 
SGXBs I Kuulkers"2005^ and the Supergi ant Fast X- ray Transients 
(SFXTs;|sguera et al. 2005, 2006, Neguer uela e7ai][2 006). 

Around 9 SFXTs are known to date. They spend most of the 
time in a low level of X-ray activity characterized by X-ray lumi- 
nosity values in the range IO^--IO'''* erg s"', well below the per- 
sistent bright state of other classical SGXBs (~I0''^ erg s"'). Only 
occasionally do they display bright and fast X-ray outbursts with 
a duration no longer than a few days (typically only a few hours) 
and a dynamical range in flux of lO^^-IO''. A number of hypothe- 
ses have been proposed to explain the peculiar flaring behaviour 
of SFXTs, both as a class of objects on their own, and in the con- 
text of persistently accreting supergiant X-ray binaries. Some of 
the models proposed invoke structure in the wind of the supergiant 
companion. This structure could be either in the form of clump- 
ing of spherically symmetric outflow from the supergiant dono r 
din 't Zandl2003 : lLevder et^bOOTl : IWalter & Zurita Heraj|2007h . 
or in the form of an equatorially density enhanced wind from the 
supergiant, incline d at some angle to the orbit of the neutron star 
dSidoli et alj2007l) . Futhermore, variations in the orbital eccentric- 



ity jNegueruela et al.ll2008l) can be used to explain the difference 
between supergiant X-ray binary systems showing transient as op- 
posed to persisten t X-ray emission. A different model proposed by 
IBozzo. Falanga, & Stella 1 2008b) makes use of transitions between 
accretion gating mechanisms, such as centrifugal and magnetic bar- 
riers, brought about by variations in the stellar wind, to explain the 
large dynamic range in flux observed in these systems. 



Fast outbursts from SFXTs are relatively rare occurances (i.e. 
a few are detected per year), and long periods of inactivity are 
interspersed here and there with short flares spaced by irregu- 
lar intervals of time. Because of this, the identification of peri- 
odicity is a very challenging task which can be fullfilled only 
through long-term and continous monitoring. To date, while re- 
current (but non-periodic) outbursts have been observed in many 
SFXTs, only one, IGR JI 1215-5952 has been shown to exhibit pe- 
riodic flarin g activity se parated by reg ular interva ls of ~ 165 days 
l lSidoh. Paizis, & Mereghettil200dlRomano et alj2007 ). The iden- 
tification of periodicity in the outburst behaviour of SFXTs, most 
likely represents the orbital period of the binary system, and is a key 
diagnostic for studying the geometry of the system (i.e. orbital ra- 
dius, eccentricity) and hence for understanding the physical reasons 
behind their very unusual X-ray behaviour. In this letter we report 
a comprehensive temporal study of the SFXT SAX J1818. 6-1703, 
leading to the discovery of a ~ 30 day periodicity, making this the 
second confirmed periodic SFXT. 
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2 PREVIOUS OBSERVATIONS 

SAX J1818.6-1703 was reported for the first time by Bep- 
poSAX/WFC on March 11, 1998 during one o f the Galactic Bulge 
observation campaigns jin't Zand etalll 19981) . The peak flux was 
~400 mCrab (9-25 keV) and the outburst lasted for 0.1 days and 
was characterised by a fast decay. 

On September 9, 2003, the IBIS telescope on INTEGRAL de- 
tected 2 short and intense outbursts with fluxes reaching 400 mCrab 
in the 18-45 keV en ergy band. This source activity was reported 
by iGrebenev & Sunv aev (2005) as showing a complex time pro- 
file lasting for about 1 day, while the main event duration was 
2.7 hours. Spectra at different evolutionary phases extending up to 
~100 keV were fitted with thermal bremsstrahlung models. The 
nature of the source was still unclear but RXTE and INTEGRAL 
data exhibiti ng such outbursts w ere considered indicative of a pos- 
sible SFXT jSguera et al.ll2005l) . a picture later strengthened by 
iNegueruela & SmitH ( l20o5 ) on the basis of an optical/NlR coun- 
terpart with an 1-band spectrum typical of an OB supergiant. Thi s 
result was confirmed soon afterwards by lin 't Zand et alj j2006t) . 
with a Chandra observation revealing only one bright and vari- 
able source within the error region of SAX J1818. 6-1703, and 
0.21" and 0.47" from the 2MASS and USNO-Bl.O positions of 
the proposed counterpart. The authors reported a strongly absorbed 
{Nh ~ 6.0 X 10"^cm"-) power law spectrum with a photon in- 
dex of 1.9 and a 0.5-10 keV flux of 7.5 x 10"'- erg cm^^ s"' and 
also suggested the source may not be a true transient but a per- 
sistent HMXB with luminous flares such as 4U 1907+097 or Vela 
X-1. Masetti et al. (2008) confirmed the OB supergiant nature of 
the counterpart and inferred a distance of 2.5kpc. 

The source was not detected by XMM during a 13 ks obser- 
vation on October 8, 2006. F ollowing the new report of activity by 
iGrebenev & SunvaevI jlOO^, a Swift ToO was activated on April 
18, 2008 and the flux value recorded was consistent with the Chan- 
dra observation, while the XMM upper limit is a factor of 80 be- 
low these detections. With these new findings. IBozzo et"al] ( l2008al) 
concluded that SAX J1818. 6-1703 shows an X-ray flux dynamic 
range similar to that observed for other SFXTs. 



3 DATASET AND ANALYSIS 

IBIS/ISGRI data is organised and analysed in short pointings (sci- 
ence windows, sew) of typically 2000s duration. We analysed 
IBlS/lSGRl data for all available science windows covering the 
date range from March 12, 2003 to March 12, 2008 (MJD range 
52710-54537) in the 18-60 keV energy band using the standard 
OSA 7 analysis tools to provide intensity and significance images 
for each science window. We excluded data flagged as taken dur- 
ing bad time intervals, and additionally removed science windows 
with high image rms. Thereafter, we extracted fluxes from each 
science window at the position of SAX J1818. 6-1703 to create a 
light curve at science window resolution. For periodicity analysis, 
science windows where SAX J1818. 6-1703 was more than 12 de- 
grees off-axis were excluded to avoid any effects due to incomplete 
off-axis response correction dSguera et alj|2007h . After filters were 
applied, the IBIS dataset consisted of 3756 sew for a total effective 
exposure of 6. 1 Ms. 

The Swift/BAT 15-50 keV fight curve for SAX J1818.6-1703 
was obtained from the BAT Transient Monitor programme 
dKrimm et al.l2006t) . and covered the date range from May 10, 2005 
to August 23, 2008 (MJD range 53500-54705). Data were filtered 
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Figure 1. Lomb-Scargfe periodogram for the IBIS/ISGRI 18-60 keV (up- 
per) and Swift/BAT 15-50 keV (lower) light curves. The 99.9% significance 
power is 15. 1 for IBIS and 14.8 for S wift/BAT. derived using the standard 
formulation iHome & Baliunaslll98d) 



to include only pointings with data quality flag (good quality). 
After filtering, 10975 individual pointings with a total effective ex- 
posure of 5.6 Ms were retained. 



4 PERIODICITY ANALYSIS 

A standard test for p eriodicity was carried out using the Lomb- 
Scargle method dLomb 1976; Scargle 1982, 1989) on the IBIS com- 
bined light curves in the 18-60 keV band. An indication of periodic 
signal at around 30 days was found, but with a significant degree of 
aliasing and a strong first harmonic (Figure[T] upper panel). 

Some aliasing is not unusual in IBIS light curves, because of 
the window function imposed by the approximately 6-monthly ob- 
servation pattern for sources near the Galactic Bulge. Furthermore, 
significant leakage of power into the harmonic is usually consid- 
ered an indicator of non-sinusoidal modulation, which the Lomb- 
Scargle algorithm must compensate for by the addition of higher 
frequency components. To provide an independent verification of 
the suspected periodicity, the same analysis was performed on the 
Swift/BAT 15-50 keV light curve (Figure [T] lower panel). In this 
case, the aliasing is much reduced due to the more even temporal 
coverage provided by the large BAT field of view, but the power of 
the periodicity is substantially reduced due to the lower sensitivity 
per pointing compared to IBIS/ISGRI. Nevertheless, the periodic- 
ity at ~30 days is still detected at above 99.9% significance in both 
light curves. 

In order to investigate further, the detected periodicity was ver- 
ified by the use of an epoch-folding method on the IBIS/ISGRI light 
curve, this method being less dependent on the periodic flux profile. 
This search again resulted in a clear maximum at around 30 days 
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Figure 2. Epoch folding analysis of the IBIS/ISGRI 18-60 keV lightcurve. 




Figure 3. IBIS light curve folded on the 30 day period. The narrow line 
is for the whole dataset, while the bold line is the result when the largest 
outburst is removed. The arrows refer to the X-ray measurements discussed 
in Sectionl?] 



(Figure [2J. Fitting the peak yields an estimate of both the period 
and its associated error, specifically 30.0 ±0.1 (Icr) days. 

The IBIS light curve folded on the 30 day period (Figure [3} 
shows a very significant modulation, and also shows that the source 
has a duty cycle of around 4-6 days during each 30 day period. 
Simulations performed on data with the same duty cycle verified 
that the strong harmonics observed in the Lomb-Scargle analysis 
are due to this profile. 

In order to ensure that the large outburst at around MJD 52891 
was not dominating the overall folded light curve, the analysis was 
repeated removing the largest outburst (specifically data for MJDs 
52889-52893). Figure [3] also shows the folded IBIS light curve 
with the outburst removed, and there is no substantial difference 
in shape with or without the large outburst. Indeed, the IBIS light 
curve analysed represents some 70 cycles of the ~30d period, and 
should be relatively robust against any such dominance by a single 
outburst. Furthermore, the analysis in section [6] shows that SAX 
J1818. 6-1703 has a high degree of recurrence, indicating that this 
folded light curve should be a reasonable representation of the av- 
erage outburst. 




53670 53671 
WD 



53851 
MJD 



Figure 4. IBIS 18-60 keV light curves showing the previously unreported 
flares at MJDs 53671 and 53851 (flares B and C in Table[T). 



5 SEARCH FOR NEW OUTBURSTS 

Knowledge of the period and phasing of the known outbursts is of 
great benefit when examining the light curves for events of lower 
significance. A search was carried out for other outbursts with sim- 
ilar characteristics to the known ones in both the IBIS and BAT 
light curves. A significance was determined for a 4-day period 
around each of the expected outburst times where sufficient data 
were available. Table [T] reports detections above 7 sigma in either 
instrument, and the corresponding significance in the other instru- 
ment lightcurve. 

The IBIS 18-60 keV light curves of the newly discovered out- 
bursts around MJDs 53671 and 53851 (labelled B and C in Table[3 
are shown in Figure|4] 

The period of flaring activity around MJD 53671 is similar to 
that observed in other SFXTs, a period of activity of around 2 days 
made up of several faster flares lasting a few hours. By contrast, the 
outburst detected around MJD 53851 shows a much cleaner struc- 
ture, with just two detectable flares. Note that, apart from the differ- 
ent intensities of the outbursts (around 50 mCrab and 120 mCrab 
respectively) the data quality for the earlier flare is also degraded 
because the source is further from the telescope axis. 



6 OUTBURST RECURRENCE 

While only detections above 7 sigma are reported in Table[T] a num- 
ber of other detections in the significance range 3-7 sigma were in- 
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Table 1. New outbursts from SAX J1818. 6-1703 detected using the period determination 





MJD 


IBIS significance 


IBIS pointings 


BAT significance 


BAT pointings 






(18-60 keV) 




(15-50 keV) 




A 


53279.8 - 53283.8 


7.5 


97 


no data 


no data 


B 


53669.4 - 53673.4 


14.8 


68 


0.5 


27 


C 


53849.3 - 53853.3 


14.5 


37 


6.5 


35 


D 


54388.8 - 54392.8 


3.9 


5 


9.9 


42 



LJ 



4-day significance 

Figure 5. Significance distributions for all available 4-day periods around 
periastron (red) and away from periastron (blue, scaled) 



dicated during the search for new outbursts. This suggests that the 
recurrence rate for outbursts is quite high for SAX J1818. 6-1703, 
and to investigate this further, the IBIS light curve was analysed 
to determine the 4-day significances both around, and well away 
from, the predicted outburst times. 

For the IBIS light curve, data were available for 23 possible 
outburst times where more than two sew fell within the 4 day win- 
dow, while 68 independent tests were performed for 4-day win- 
dows whose starts were offset by between 9 and 17 days after a 
predicted outburst. Of the 23 tests performed around the predicted 
outbursts times, 15 produced a significance above 3 sigma, while 
none of the 68 tests away from those times did. The distribution 
of 4-day IBIS significances is shown in Figure [5] and it is evident 
that a large fraction of the predicted outbursts produce emission 
above the quiescence level. Furthermore, at least two of the appar- 
ent non-detections may be due to limited exposure during potential 
outbursts, and so a conservative lower estimate can be made that 
detectable emission is seen on more than 50% of the times when an 
outburst is predicted by the 30-day period. 



7 DISCUSSION 

This 30 day periodicity, evidenced by a regular increase in emission 
and occasional strong outbursts, is strongly indicative of a binary 
orbital period. The various observations of SAX J1818. 6-1703 per- 
formed to date, placed in the context of this new orbital period, 
are summarised in Table |2] This table includes both the INTE- 
GRAL/IBIS and Swift/BAT detections, together with the X-ray ob- 
servations performed by various telescopes. 

The results of the X-ray observations carried out to date are 



consistent with this orbital period, and for reference the phase (ar- 
bitrarily referenced to the BAT trigger at MJD 54540.659) of the 
various X-ray observations are indicated in Figure[3l It can be seen 
that the XMM observation, which reported a firm upper limit on the 
flux, occurred at a phase of 0.51, i.e. apastron. Thus this measure- 
ment can be interpreted as setting an upper limit on the quiescent 
flux from the system. By contrast, the Chandra observation was 
performed at phase 0.89, and the reported monotonic rising flux 
appears to show the onset of emission, but whether this is a slow 
increase as the compact object approached periastron, or the onset 
of an outburst, we cannot say. The XRT observation was performed 
at a phase of 0. 1 1 and indicates some low-level activity, compatible 
with the Chandra flux, and an order of magnitude above the XMM 
upper limit. 

We performed a similar analysis using the RXTE/ASM light 
curve. Unfortunately, while we can conclude that RXTE/ASM is 
detecting X-ray emission from the source, it does not provide any 
further diagnostic information due to its limited instantaneous sen- 
sitivity to a source that is only rarely above 20mCrab, and time 
coverage that does not really manage to capture the detail of the 
flaring time structure. 

A systematic search of the IBIS data for both the large out- 
burst and newly discovered outbursts has been performed to look 
for pulsations in the range 0.2 - 5000s. No statistically significant 
pulse period has been detected by this analysis and indeed no pulsa- 
tions from SAX J1818. 6-1703 have been reported in the literature. 

We conclude that although SAX J1818.6- 1703 was originally 
detected as a transient, and later categorised as a SFXT, knowledge 
of the orbital solution enables us to see a rather different picture. 
SAX J1818.6-1703 is a supergiant system with a compact object 
in a 30 day orbit, and shows regular emission on at least 50% of the 
periastron passages. Moreover, some periastron passages give rise 
to larger outbursts, during which fast flaring structure is evident in 
the light curve - it is this structure which has led to the identifi- 
cation as a SFXT. SAX J1818. 6-1703 shows an extreme dynamic 
range of emission, from < 1 . 1 x 10"" erg.cm"-.s"' (0.5-10 keV) at 
apastron jBozzo et al.l2008al) . to 3.7x10"' erg.cm"^.s"' by extrapo- 
lation of the best-fit spectrum during the largest recorded IBIS out- 
burst ( Grebenev & Sunyaev 2005) to 0.5-lOkeV, which again sup- 
ports the SFXT categorisation. The simplest interpretation of this 
behaviour is of a compact object in an eccentric orbit, periodically 
accreting from the wind of a donor supergiant. 

As the second SFXT to be shown to have a periodic out- 
bursting behaviour, SAX J1818. 6-1703 adds to our overall pic- 
ture of the supergiant binary systems, placed between the persis- 
tent systems with circular orbits, and the long-period system IGR 
Jl 1215-5952 which shows outbursts on a 165 day period. This in- 
dicates that orbital characteristics may play an important role in 
the accretion mechanisms of these systems, much as they do in 
Be transient systems. The shorter period of SAX J1818. 6-1703 
has enabled us to determine that the system shows a high degree 



30 day period in SAX J1818.6-1 703 5 

Table 2. Timeline of X-ray and liard X-ray observations of SAX J1818.6-1703 



Instrument 


MJD 


Phase 


IBIS peak flux 
(18-60 keV, mCrab) 


Notes 


References 


SAXAVFC 


50883.86 


0.07 




Discovery 


ijL't Zandetal. (1998^ 


IBIS 


52891.6 


0.01 


400 


Large outburst 


Grebenev & Sunvaev (200^ 


IBIS 


52922.96 


0.06 


50 


Flaring activity 


Seuera et al. (2005) 


IBIS 


53281 


1.00 


<50 


New 




IBIS 


53671 


0.00 


50 


New 




IBIS -1- BAT 


53851 


1.00 


120 


New 




Chandra 


53997.5 


0.89 




Rising flux 


in 't Zand et al. (2006^ 


XMM 


54016 


0.51 




Upper limit 


Bozzo et al. (2008a) 


IBIS + BAT 


54390 


0.98 


<50 


New 




BAT 


54540.659 


1.00 




BAT trigger 


Barthelmv et al. (2008) 


IBIS 


54568 


0.91 




Rising flux 


Grebenev & Sunvaev (2008) 


Swift/XRT 


54574 


0.11 




Detection 


Bozzo et al. ('2008a) 



of recurrence. The system appears similar to IGR J18483-0311, 
w hich has an orbital period of 18.52 days, spin period of 21.0526 
s iSgueraetal.ll2007l), and has recentl y been confirmed as a su- 
pergiant system ( iRahoui & Chat\l2008h . Other SFXT systems may 
eventually reveal similar behaviour, but with more eccentric orbits 
and/or lower recurrence rates which have so far prevented the de- 
termination of their orbital periods. 

We encourage detailed follow-up observations based on 
this orbital period, such as have proved so effective for IGR 
Jl 1215-5952. A detailed picture of the temporal structures seen 
during both faint and bright periastron passages should prove an 
invaluable probe into the environment from which the compact ob- 
ject is accreting. 
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